Regular Black Hole in General Relativity Coupled to Nonlinear Electrodynamics 
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The first regular exact black hole solution in General Relativity is presented. The source is a 
nonlinear electrodynamic field satisfying the weak energy condition, which in the limit of weak field 
becomes the Maxwell field. The solution corresponds to a charged black hole with \q\ < 2scm « 
0.6 m, having the metric, the curvature invariants, and the electric field regular everywhere. 
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In General Relativity the existence of singularities appears to be a property inherent to most of the physically 
relevant solutions of Einstein equations, in particular, to all known up-to-date black hole exact solutions Q. The 
Penrose cosmic censorship conjecture states that these singularities must be dressed by event horizons; no causal 
connection could exist between the interior of a black hole with the exterior fields, thus pathologies occurring at the 
singular region would have no influence on the exterior region, and the Physics outside would be well behaved (cf. |^ 
for a review on the recent status of this conjecture). 

To avoid the black hole singularity problem, some regular models has been proposed . All of them have been 
referred to as "Bardeen black holes" Q , since Bardeen was the first author producing a surprising regular black hole 
model 1^. No one of these models is an exact solution to Einstein equations; there are no known physical sources 
associated with any of them. The attempts to solve this problem have usually been addressed to the search of more 
general gravity theories. The best candidate today to produce singularity-free solutions, even at the classical level. 



due to its intrinsic non-locality, is string theory |10|. There are examples in other contexts, for instance, in = 1 
supergravity domain wall solutions with horizons but no singularities have been found (cf . pH | , and references therein) , 
another example is given in exact conformal field theory . 

We show in this Letter that in the framework of the standard General Relativity one can find singularity-free 
solutions of the Einstein field equations coupled to a suitable nonlinear electrodynamics, which in the weak field 
approximation becomes the usual linear Maxwell theory. Previous efforts on this direction with nonlinear electrody- 
namics either have been totally unsuccessful or only partially solve the considered singularity problem p3|-p^. We 
propose a new nonlinear electrodynamics which coupled to gravity actually produces a non-singular exact black hole 
solution satisfying the weak energy condition. 



The gravitational field of our solution is described by the metric 
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while the associated electric field E is given by 
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Notice that this solution asymptotically behaves as the Reissner-Nordstrom solution, i.e., 

-gtt = 1 - 2TO/r + q^jr^ + ©(l/r^), E = g/r^ + 0(l/r3), 

thus the parameters to and q are related correspondingly with the mass and the electric charge. For a certain range 
of the mass and charge our metric (|^) is a black hole, which in addition is regular everywhere. Accomplishing the 
substitutions x = f/l?! and s — |(j|/2to, we rewrite gu as 

- gu = .) ^ 1 - + (Y^^, (3) 
which, for any nonvanishing value of s, has a single minimum; cf. Fig. ^ There exists a single real critical value of x, 
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FIG. 1. Behavior of —gtt for different values of cfiarge. 



Xc, and one of s, Scj to be determined from A{xc, Sc) — and dxA{xc, Sc) = 0, namely 

-f + - - 1 = 0, 2t^ +4 = 0, 

s s s s 

where = a;^ + 1. To solve these equations, one substitutes s = t{P — 3)/(2i^ — 4) from the second equation into 
the first one arriving at — 4i* + 2t^ — 1 = 0, which has only one real solution for i^, thus the corresponding critical 
values are Sc ~ 0.317 and Xc ~ 1.58. For s < Sc the quoted minimum is negative, for s — Sc the minimum vanishes, 
and for s > Sc the minimum is positive. Evaluating the curvature invariants R, Rf^^R^'^ , and R^uapR^'^"'^ for metric 
(|l|) one establishes that they are all regular everywhere, cf. Fig. ||; thus for s < Sc the singularities appearing in (|l|) 
due to the vanishing of A are only coordinate singularities describing the existence of horizons, consequently, we are 
in the presence of black hole solutions for \q\ < 2scm « 0.6 m. For these values of mass and charge we have, under 




FIG. 2. Regular behavior of the Ricci, q^R, Ricci square, q^R^^R^" , and the Riemann square, q"^ R^.^apR'^'^"^ , scalars for 
different values of charge; the abscissa is r/\q\. 



the strict inequality \q\ < 2sc77i, inner and event horizons for the Killing field k = d/dt, defined by the real solutions 
of the quartic equation —ki^i_k^ = A = 0, which are given by 
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g{s) = 4 (9 s + 74 + \/27(400 - 112 + 47 s2 - 4) 



For |q| = 2scm, the horizons shrink into a single one, corresponding to an extreme black hole (Viy(fcpfc^) = 0). 
The extension of the metric beyond the horizons r± becomes apparent by passing to the standard advanced and 
retarded Eddington-Finkelstein coordinates, in terms of which the metric is smooth everywhere, even in the extreme 
case. Following step by step the procedure presented in ||l^, Chap.V] to derive the global structure of the Reissner- 
Nordstrom black hole, one can arrive at the global structure of our solution and construct the Penrose diagrams; 
nevertheless, because of journal length restrictions we omit here the corresponding calculations and diagrams, leaving 
this issue for an extended publication. Briefly, what one encounters in the case of our non-extreme black hole solution, 
|g| < 2scm, is the splitting of the space-time into three regions, I: r > r+, II: r_ < r < r+, and III: < r < r_; 
cf. Fig. ^ In each region one introduces advanced and retarded coordinates u and v, related with r through the so 
called tortoise coordinate r* = J A^^dr, which in our case is quite involved. Further, by the inversion of u and d, 
u — > — u, V — !■ —V, one obtains the remaining regions I , II , and III . Introducing a new set of null coordinates one 
arrives at the maximal extension of the non-extreme black hole. The Penrose diagram of the maximal analytical 
extension of our solution is obtained by gluing appropriately copies of these six regions upward and downward ad 
infinitum. In the extreme black hole case, |g| = 2scm, there arise two regions, I: r > Tc and III: < r < re, cf. Fig. 
in which again one introduces advanced and retarded u and v coordinates to accomplish the maximal analytical 
extension; these two region determine the main building block of the extension. To construct the Penrose diagram of 
the maximal analytical extension, one glues copies of this block in a suitable way. In both cases, extreme and non- 
extreme, there is no singularity at r = 0, which is now simply the origin of the spherical coordinates. Summarizing, 
our space-time possesses the same global structure as the Reissner-Nordstrom black hole except that the singularity, 
at r = 0, of this last solution has been smoothed out. 

For |g| > 2scTO, there are no horizons and the corresponding exact solution represents a globally regular space-time. 
It is worthwhile to mention in this respect the existence of globally smooth solutions to the Einstein+matter (Yang- 
Mills, Yang-Mills-Higgs) equations; although there are demonstrations of the existence of these solutions |jl^,Q, they 
are numerically given and there are no analytical closed expressions for them p^ ; cf. ]2C[ |, and references therein. 

The fields (0) and (|) arise as a solution of the Einstein-nonlinear electrodynamic field equations derived from the 
action proposed in Einstein-dual nonlinear electrodynamic theory [ pl| , which in the studied case becomes 

S^l^v[±.R-l-C(F)). (5) 

where R is scalar curvature, and £ is a function of F = jF^^F^^'^ . Alternatively, one can describe the considered 
system using another function obtained by means of a Lcgendre transformation |pl| : 

n = 2FCF - C. (6) 

Defining P^^ = CpF/^^^, it can be shown that 7i is a function of P = jP^iuP^'^ ~ {^f)^F, i.e., dH ~ 
{CF)~^d{{CF)'^F) — TipdP. With the help of Ti, one expresses the nonlinear electromagnetic Lagrangian in the 
action (|5|) as £ = 2P'Hp — TC, depending on the anti-symmetric tensor P^jy. The specific function TC, determining the 
nonlinear electrodynamic source used, is given as 



where s — \q\/2m and the invariant P is a negative quantity. The corresponding Lagrangian occurs to be 

(l-8V-2g2p-6g2p) 3 {^2q^P)^/'^ - 2^/ -2q^P ^ 
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The function (Q) satisfies the plausible conditions, needed for a nonlinear electromagnetic model, of (i) correspondence 
to Maxwell theory, i.e., Ti P for weak fields (P <C 1), and (ii) the weak energy condition, which requires 7i < 
and Hp > 0; cf. Fig. |[ We would like to point out that our solution, in addition to being regular and to satisfying 
the weak energy condition, is characterized by another feature: it does not admit a Cauchy surface. Hence, it 
does not contradict the Penrose singularity theorem supported on the hypotheses of: fulfillment of the null energy 
condition, existence of a noncompact Cauchy surface, and existence of a closed trapped surface and concluding no 
null geodesically completeness of the space-time. 
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FIG. 3. Behavior of q^Ti. and Tip with respect to the positive abscissa —2q^P for different values of charge. 

In what follows we shall briefly give the main lines of the integration process yielding the studied solution. The 
Einstein and nonlinear electrodynamic equations arising from action (|^) are 

g;: = 2{npP,,xP''^ - 5;{2Pnp - n)), (9) 



V^P"^ = 0. (10) 
In order to obtain the solution (Q), (|^), wc consider the static and spherically symmetric configuration 

5 - - - - + ^ V*^ + fi - V + ^) " '^^^ + (11) 

as 



and the following ansatz for the antisymmetric field P^j, — 2(5? (5^,1? (r). With these choices the equations (|1C|) integrate 



where we have chosen the integration constant as q since, as it was previously anticipated, it actually plays the role 
of the electric charge. The evaluation of the electric field E = Ftr — HpD, using expression (^ for H, gives just the 
formula (||). The j* component of Einstein equations (||) yields the basic equation 

"-^^ = mp). (13) 

Substituting Ti from (|^) with P = —q'^/2r'^ one can write the integral of (13) as 
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the integrand above can be expressed as dy(2my^ / {y^ + g^)'^/^ — if" /{ip' + q^Y), thus one arrives at 

n 4 2 4 

n o 2mr* q 

g = 2mr- p-j— 5^ + p-^— (15) 
Substituting Q into —gtt = 1 — 2m /r + Q/r^ one finaUy gets Eq. (|l|). 
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